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The possibihty of using single molecule junctions as electron pumps for energy conversion and 
storage is considered. It is argued that the small dimensions of these systems enable to make use 
of unique intra-molecular quantum coherences in order to pump electrons between two leads and 
to overcome relaxation processes which tend to suppress the pumping efficiency. In particular, we 
demonstrate that a selective transient excitation of one chromophore in a bi-chromophoric donor- 
bridge-acceptor molecular junction model yields currents which transfer charge (electron and holes) 
unevenly to the two leads in the absence of a bias potential. The utility of this mechanism for charge 
pumping in steady state conditions is proposed. 



Single molecule junctions were explored and studied 
extensively in recent years[l] primarily due to their po- 
tential for nano-electronics devices. Indeed, non triv- 
ial current-voltage characteristics were measured, reflect- 
ing the unique properties of organic molecules, associ- 
ated with their detailed structure [2-4]. Less attention 
was given to the possibility of using single molecule 
junction as energy conversion devices. However, the 
successful implementation of organic molecules in solar 
cells [5], and the quest for energy conversion units on 
the nanoscale, raise an intriguing challenge of design- 
ing particular molecular structures that can function as 
energy conversion units in a junction architecture. In 
this work bi-chromophoric molecules are suggested as 
appropriate candidates for this application. Asymmet- 
ric Donor-Bridge- Acceptor molecules in which two chro- 
mophores are separated by a molecular bridge were sug- 
gested in the past as effective molecular current rectifiers 
[6]. Below we study their function as electron pumps in 
a single molecule junction. The small dimensions of the 
proposed systems enables to make use of their unique 
intra-molecular quantum coherences in order to pump 
electrons between the leads and to over come relaxation 
processes which tend to suppress the pumping efficiency. 
Current generation by electron pumping is facilitated in 
the absence of bias by a time-dependent excitation. The 
possibility of controlling electronic currents in molecu- 
lar junctions by periodic external fields was discussed 
intensively[4, 7-10]. In particular, monochromatic, bi- 
chromatic[4, 10] or pulsed [11-13] electromagnetic radi- 
ation fields were proposed as means of controlling the 
current intensity under steady state conditions. Given 
an appropriate asymmetry of the junction such exter- 
nal fields can induce net current generation even in the 
absence of bias. This was anticipated, e.g, when the cou- 
pling to the electrodes differs in the ground and excited 
molecular state [14]. Although pumps are usually dis- 
cussed in the context of periodic driving [15-17], in order 
to pinpoint the underlying mechanism of the molecular 
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energy converter, we consider below a single pumping cy- 
cle in which one of the chromophores (termed " donor" ) 
is subjected to a sudden excitation. 

The DBA molecule is represented as three electronic 
sites corresponding to the donor(D), bridge(B) and ac- 
ceptor(A) molecular groups (see Fig. 1). Invoking for 
simplicity a non interacting spinless electrons picture, the 
molecular model Hamiltonian reads. 
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dl^ i{djn,i) is a creation (annihilation) operator for an 
electron at the / th single particle state (orbital) of the 
m th molecular site, Sm.i is the corresponding orbital 
energy, and {/3^ ; } are the inter-site coupling (hop- 
ping) integrals between the two manifolds of orbitals at 
the m th and m-1 th neighboring sites. The indexes 
I = 1 and I = 2 point to the Highest Occupied Molecular 
Orbital(HOMO) and Lowest Unoccupied Molecular Or- 
bital(LUMO) at each site, and the indexes m — 1,2,3 
represent the D, B and A sites respectively. In view 
of the relatively large HOMO-LUMO energy gaps, the 
inter-site couplings were restricted to a diagonal form. 
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FIG. 1. A scheme for DBA molecular junction at zero bias. 
Notice that the LUMO levels of the D and A sites are degen- 
erate but the HOMO-LUMO energy gaps are different. An 
initial D excitation (marked by the dashed arrow) is followed 
by transient electronic currents (thick arrows) from left to 
right. 
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l^mi i' ~ Pm,i5i i> (off-diagonal terms yielded insignificant 
changes to the results reported below), and all inter-site 
couplings were given the same value, /3„i_/ = (3. The 
value of /3 is strongly sensitive to the bridge length and 
composition, and determines the time scale for coherent 
tunneling between D and A. Since coherent tunneling 
plays an important role here, the value was taken larger 
than KbT. Asymmetry of the bichromophoric molecule 
is reflected in different HOMO-LUMO gaps at the donor 
and acceptor sites, i.e., £i_2 — £1,1 > £3,2 — £3,1, where 
the effective electronic coupling between the D and A 
sites is mediated by the bridge HOMO-LUMO gap. Here 
we focus on a particular system in which the D and A 
LUMO states are degenerate, i.e £12 = £3.2- (a similar 
effect would be obtained if the degeneracy would apply 
to the HOMO states). A schematic representation of the 
model is given in Fig.l. In view of the differences between 
the HOMO-LUMO gaps at the different molecular sites, 
direct intra-molecular exciton energy transfer [18] is in- 
efficient, and terms of the type (oc d}^^ i''^m' i) 
were excluded in Eq.l. 

The terminal sites (D and A) are coupled to two (left 
and right, respectively) electrodes, so that the full Hamil- 
tonian reads. 
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The electrodes are modeled as reservoirs of non inter- 



acting electrons, He 
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where 



'^fc ji^kj) is a creation (annihilation) operator of an elec- 
tron in the fcth single particle state of the Jth elec- 
trode' conduction band. The molecule-electrodes cou- 
pling terms correspond to electron transfer (hopping) 
between the D and A sites and the respective elec- 
trode, i.e., He,M = Vl ® Ul + Vr ® Ur + h.c, with 

and Ur = J2k '^k,Rbk.R- For simplicity, electron en- 
ergy transfer from the molecular excitons to electron- 
hole pairs creation at the leads [19] is not accounted for. 
The values of {uk,j} reflect the spectral density func- 
tion of the Jth electrode. Below, each electrode is mod- 
eled as a half-filled tight binding chain with a bandwidth 
|47j| and an electrochemical potential, fjij. The elec- 
trode band is discretized [20], where the single particle 
energies are Sk.j = hj — 2\^j\cos[-^^], and therefore, 

Mfc,j = £.J\J Wyi^^'^^WTi^- hopping parameter 

between the terminal molecule and electrode sites. 

An initial sudden preparation of an exciton at the 
molecular donor site is represented by a factorized den- 
sity operator, p{0) = pm{0) ® pr ® pl ■ Pm{^) is the 
molecular density which accounts for the population of 
electrons {n^.i = 1) or holes (rtm,i = 0) in the molecular 
single particle states (see Fig.l and Table I), pm{0) ~ 
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TABLE I. Model Parameters (Energy values are in eV) 
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The reservoirs 
with equilibrium 

For weak coupling between the molecular system and 
the electrodes, it is convenient to regards the reservoirs 
as external baths, and follow the system time-evolution 
using a Redfield type equation of motion for a reduced 
system density [21], pMit) = ^''[/o(^)]j keeping terms up 
to second order in idg jv,/[20], 

d i • 
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electrode correlation functions, 

EJ"fc,j|'e-^^^'''"[l-/,/(£fe.j)] and 

Ek M'ei^^-'^fjiskj). fAsk.j) ^ ^ 
is the electrode Fermi function. 



The change in the charge on the molecular sys- 
tem is defined as, = e • tr[N^pM{t)], where 



dt y dt^ 

= Y,m=i SLi '^ln,i'^rn,i is the electronic number op- 
erator in the system. Using [N, Hm] = 0, one can express 
the charge dynamics in terms of additive contributions 
from the two electrodes, i.e., = lR{t) + lL{t), 

where the transient current from the electrode into the 
molecule read, 

Ijit) = AeRe{tr{N[{Fj{t)pM{t) - /5M(t)F](t)}, Vj]}) (5) 

The integral, Qj — Ij{t)dt, is the accumulated 
charge due to electron and hole transfer (positive val- 
ues correspond to electron transfer into the molecule), 
which depend on the intramolecular dynamics, including 
intramolecular coherence. For the DBA molecule illus- 
trated in Fig.l, the HOMO energies at the D and A sites 
are detuned, such that the effective electronic coupling 
between them is week regardless of the B HOMO energy. 
In contrast, the two LUMO energies are resonant and 
therefore the electronic coupling between them depends 
sensitively on the B LUMO energy. 
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FIG. 2. (a) Accumulated charges at the electrodes (current 
integrals) as functions of the bridge LUMO energy. Inset: 
the transient currents for a selected energy, £2,2 = 0.05eV. 
(b) Transient populations of the molecular states following 
a sudden excitation of the donor for £2,2 ~ 0-05eV. Each 
molecular orbital is denoted by its site (D, B or A) and energy 
indexes (1 for HOMO and 2 for LUMO). 



In Fig. 2 the accumulated charges are plotted as func- 
tions of the bridge LUMO energy, £2,2 [22] (the model 
parameters are detailed in Table I ). As one can see, 
there is a net charge transfer from the left electrode to 
the right one, > 0, Qr < 0, (see also the inset, for the 
transient currents). The "charge pumping" efficiency, re- 
garded here as the net amount of charge transferred from 
the left electrode to the right one, increases and reaches a 
maximum as the B level approaches resonance with the 
D and A levels. The charge transfer between the elec- 
trodes reflects intramolecular dynamics as presented in 
Fig. 2b. The apparent dominant processes following the 
initial excitation at the donor are a decay of the hole 
at the donor HOMO (see D,l), and two decay processes 
of the electron from the donor LUMO: One is a direct 
relaxation to the left electrode (see D,2) and the other 
is an indirect relaxation to the right electrode following 
coherent tunneling from the D to the A LUMO level (see 
A,2). When the bridge LUMO is detuned from the D/A 
LUMOs (i.e., when |/3| << |e2,2 — £1,2!) the period of the 
DA tunneling oscillations can be approximated by the 
McConnell formula [23], tda « ^^^^^y^jP^, and becomes 
shorter as the gap between the B and the D(A) orbital 
energies decreases. As the period becomes sufficiently 
short on the time scale of the electron relaxation to the 
leads, electron hopping to the right becomes as frequent 
as electron hopping to the left. Since the hole relaxes pre- 
dominantly to the left lead, the balance of these processes 
is net transfer of negative charge from left to right. 

It is emphasized that although the simple model 
outlined above is formulated in the language of non 
interacting particles, it can guide a search for real 
bichromophoric systems where the proposed mechanism 
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FIG. 3. Accumulated charges at the electrodes (current inte- 
grals) as functions of the bridge reorganization energy. 

of coherent pumping following a local excitation of a 
donor chromophore should be observed. Considering a 
bichromophoric molecule where the two different chro- 
mophores, are weakly coupled to each other and to the 
leads, the condition for near degeneracy of the two LU- 
MOs in the non-interacting model corresponds to a gen- 
eralized condition for a near degeneracy between the first 
excited (many body) state of one of the chromophores the 
" Donor" and a charge transfer (many body) state of the 
bichromophoric system. In such a case, local excitation 
of the Donor would result in to intra-molecular dynam- 
ical charge oscillations, followed by the proposed charge 
pumping. 

Let us now consider a possible effect of bridge vibra- 
tions on the pumping efficiency. For this purpose we ex- 
tend the model to account for linear on-site coupling of 
the bridge excited state to a harmonic mode. The model 
DBA Hamiltonian (Eq.(l)) is replaced by, 

HM.v^b = Hm + \ {c) + C)4 2^2,2 + hn{c^ C + 1/2) (6) 

v2 

The last term is the nuclear contribution, where (c) is 
a creation (annihilation) operator of a vibration quan- 
tum. Ml, and the vibronic coupling strength is mea- 
sured by the bridge reorganization energy, A = 2hn- 
In Fig. 3, the effect of A on the pumping efficiency is 
demonstrated. The acceptor and bridge HOMO states, 
which remain occupied and hardly participate in the 
dynamics (See Fig. 2), were excluded in this case (see 
Table I), and the initial electronic excitation was rep- 
resented as pM,vtb{^) = Pm{0) <8) |0)(0|, where (|0)) is 
the vibrational ground state. The left-to-right electron 
pumping efficiency is shown to modulate as a function 
of the vibronic coupling strength, where efficiency peaks 
{—Qr = Ql ~ 0.5) are obtained near A w 0.05,0.11. 
As above, the increased efficiency corresponds to reso- 
nances between the donor and acceptor LUMO orbitals 
and, in this case, vibronic eigenstates of the bridge, where 
the peak efficiency corresponds to £2.2 — A « 0, hfl re- 
spectively. Again, approaching the resonance condition 
guarantees that the frequency of tunneling oscillations 
between the D and A exceeds the rate of charge relax- 
ation to the electrodes, and thus the charge pumping 
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efficiency increases. Notice that the efficiency peaks in 
Fig. 3 are characterized by a slow rise, followed by a fast 
fall as A increases. This reflects a characteristic Franck- 
Condon blockade effect [3]. Focusing on the first reso- 
nance, one can approximate the tunneling frequency by 

a generalized McConnell formula [24] , v = "^^^^I^^ , (for 

£1.2 — £3,2 = 0). As A rises from zero, two effects are 
competing. The energy barrier for through-bridge tun- 
neling (|e2,2 — A|) decreases, and the hopping frequency 
to the bridge, scaled by the Franck Condon (FC) over- 
lap between the ground vibrational states of the empty 
and electronically charged bridge, decreases as well. For 
sufficiently small A the exponential drop in the FC fac- 
tor is minor, and the decrease in the barrier energy leads 
to increasing tunneling frequency. Increasing A farther, 
beyond the resonance point, the tunneling frequency de- 
creases due to a rise in \e2,i — A|, but the fall is much 
sharper than the rise since the donor-bridge hoping is 
additionally blocked as the FC factor becomes exponen- 
tially small. A similar qualitative modulation is observed 
for larger A values due to the first excited vibration state 
of the bridge. 

In conclusion, the versatility of organic compounds 
and the ability to control their electronic structure by 
chemical substitutions should make the preparation of 
bi-chromophoric DBA molecules as discussed here a re- 
alistic task. The advances in fabrication of molecular 



junctions (of single molecules or of ordered monolayers) 
suggest that devices can be based on the single molecule 
properties. In particular, irradiation of one of the chro- 
mophores can be converted into charge pumping even in 
the absence of a bias potential. The focus of the present 
work was on a single pumping cycle. An extension of this 
study to account for the details of the excitation field 
would enable to consider the operation of a multi-cycle 
(periodic) electron pump at steady state. According to 
the above analysis, the pumping direction would refiect 
the selection of either one of the two chromophores as 
the "donor" by the excitation field. Moreover, the pos- 
sibility for pumping at steady state against an opposing 
bias voltage suggests the utility of the proposed molecular 
junction for electric energy storage. The analysis of the 
appropriate field parameters that would induce the re- 
quired functions given the particular molecular electronic 
structure is a future challenge, as well as the inclusion 
of several effects that were left out for simplicity in the 
present work, but are not expected to change its qual- 
itative conclusions. These include the role of Coulomb 
repulsion and electronic correlations at the molecule, the 
vibronic structure of the chromophores, and the effect of 
energy transfer to the leads. This research was supported 
by the Israel-U.S. Bi-national Science Foundation and by 
the Israel Science Foundation. 
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